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ABSTRACT 

We compare the peculiar velocities of nearby SNe la with those predicted by 
the gravity fields of full sky galaxy catalogs. The method provides a powerful 
test of the gravitational instability paradigm and strong constraints on the 
density parameter (3 = tt°- 6 /b. For 24 SNe la within 10,000 km s _1 , we find the 
observed SNe la peculiar velocities are well modeled by the predictions derived 
from the 1.2 Jy IRAS survey and the Optical Redshift Survey (ORS). Our best 
(3 is 0.4 from IRAS, and 0.3 from the ORS, with f3 > 0.7 and (3 < 0.15 ruled 
out at 95% confidence levels from the IRAS comparison. Bootstrap resampling 
tests show these results to be robust in the mean and in its error. The precision 
of this technique will improve as additional nearby SNe la are discovered and 
monitored. 

subject headings: sup ernovae: general ; cosmology:observations, large-scale structure of 
universe — Local Group 
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1. Introduction 

The motions of galaxies probe the size of potentials formed by the gravitating matter in 
the Universe. Redshift surveys map the spatial distribution of luminous matter. Together, 
the two measure the contribution of luminous and dark matter to the contents of the 
Universe, providing a direct measure of the density parameter on the largest possible scale. 

Previous attempts to test the gravitational instability paradigm and constrain the 
mass density parameter, Q , are hindered by imprecise distance estimates to individual 
galaxies. Distance indicators based on empirical relationships between galaxy luminosity 
and internal velocity (i.e., Tully-Fisher and D n — a) yield individual distance uncertainties 
of 20-25% (Jacoby et al 1992, Tully & Fisher 1977, Dressier et al 1987, Willick et al 1996). 
To increase the utility of such data, workers have amassed it in great quantity. While this 
strategy diminishes the random component of error, imprecise distances give rise to more 
troublesome systematic errors including selection bias, Malmquist bias, and smoothing 
biases (Strauss and Willick 1995). The antidote for such biases depends on knowledge of 
the distance uncertainty, a controversial quantity (Willick 1995, Mathewson & Ford 1994, 
Federspiel, Sandage, & Tammann 1994, Giovanelli & Haynes 1996). The forced marriage 
of inhomogeneous catalogs could result in additional errors in the inferred velocity field. 
Despite these challenges, remarkable progress has been made in this field (Dekel 1994, 
Dekel, Bertschinger, & Faber 1990, Davis, Nusser, & Willick 1996, Hudson 1994, Shaya, 
Peebles, & Tully 1995, Nusser & Davis 1994). 

Type la supernovae are well suited to provide an independent test of the gravitational 
instability paradigm and to constrain the mass density. Light curves from the Calan/Tololo 
Survey (Hamuy et al 1993, 1996) and the CfA survey (Riess 1996, Riess et al 1997) 
yield distances with 5-10% uncertainty over the redshift range 1000 < cz <36000 km 
s -1 . Although the sample of observed SNe la is relatively small, the depth and precison 
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of supernova distances provide some advantages. For the reduction of random errors, 
one SN la is worth ~ 10 Tully-Fisher or Dn-a measurements. Systematic bias, which 
rises with the square of the distance uncertainty, is also ~ 10 times smaller for SNe la 
distances. Individual distance errors can be derived from SNe la light curves (Riess, Press, 
and Kirshner 1995a, 1996) which give meaningful measurements of the deviation from 
smooth Hubble flow. Combining the SNe la data with models of the predicted "peculiar 
velocities" , one can constrain the mass density traced by galaxy fluctuations while testing 
the gravitational instability model for structure formation. 



The direct comparison of observed peculiar velocities and mass density fluctuations is 
not trivial. The problem becomes tractable only by assuming that the large scale component 
of the flow is single-valued and irrotational. Since galaxy fluctuations on scales larger than 
10/i -1 Mpc are observed to be less than unity (Davis and Peebles 1983), these assumptions 
are reasonable. Furthermore, for such fluctuations at late times, linear perturbation theory 
is quite accurate, and the expected velocity field v p can be related to the peculiar gravity 



This equation simply states that the observed peculiar velocity results from the gravitational 
acceleration acting over a Hubble time. 

The gravity field g(r) may be inferred from the distribution of galaxies which are 
assumed to trace the matter field. It is common to assume an unknown but linear bias 
exists between galaxy fluctuations 5q and matter fluctuations 5 P , by employing a bias 
parameter b, i.e. 5g = b5 p . This simplification, while not accurate on small scales or in 



2. Analysis 



field g (Peebles 1980) by 
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dense environments, should suffice on large scales. Using linear theory, one can constrain 
only a combination of fi and b, (3 = VL°q % /b. A variety of methods for estimating g(r) 
from galaxy surveys are reviewed by Strauss and Willick (1995). Here, we use the method 
of Nusser and Davis (1994), which can successfully repair the gravitational distortions in 
redshift space. However, this technique cannot properly treat regions with multivalued 
relations between velocity and distance, expected near clusters of galaxies. The minimum 
smoothing scale of the derived gravity field is 500 km s _1 , and this smoothing increases 
with distance in proportion to the mean interparticle spacing of the galaxy catalog. This 
increase is necessary to suppress artificial two-body acceleration (Strauss et al 1992). 

For comparison to the SNe la peculiar velocities, we computed the gravity field from 
the 1.2 Jy flux limited IRAS redshift survey of galaxies (Fisher et al 1995). This catalog 
contains 6010 galaxies with median redshift of 6000 km s _1 and provides a useful measure of 
the density field out to cz 15, 000 km s _1 . The IRAS catalog is known to underrepresent 
the elliptical-rich cores of clusters relative to optically catalogs. Thus it is of interest to use 
the ORS catalog of optically selected galaxies for an alternative calculation of the gravity 
field (Santiago et al 1995, Baker et al 1997). The ORS contains 8457 galaxies, supplemented 
by the 1.2 Jy IRAS survey in sky regions without optical coverage. 

We measured the distances and their uncertainties for 25 SNe la within 10,000 km s _1 
using MLCS distance measurements (Riess, Press, and Kirshner 1996). The redshifts of the 
SNe la come from Hamuy et al (1996), Riess (1996), and Riess et al (1997). Because of the 
limitation of linear biasing and the risk of multivalued flows, we conservatively discarded 
one SN la, SN 1992G, due to its proximity to the Virgo Cluster. In the directions of the 
remaining SNe, we computed the expected distance-redshift relations based on the gravity 
fields of either the ORS or IRAS surveys, as a function of f3. Figure 1 shows an example of 
these curves. We assume a redshift error of 200 km s _1 as an estimate of the small scale 
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component of the radial peculiar velocity not describeable by linear theory 

For each SN, we have the distance d{ (in km s -1 ), a distance error Sdi (in km s -1 ), a 
redshift z i: and a "redshift error", 5czi = 200 km s -1 . From the gravity maps, we have the 
functions z(j,(3) and d{j,(3) along a set of points j towards each SN's direction. We seek 
the minimum separation of each point from the predicted curves, in units of the standard 
deviations. That is, for the ith SN la, we compute a contribution to a x 2 , 

r f o\2 ■ I f d i- d {h P) Y , f zj- z(j, (3) \ 2 \ 

8 X {p\ = mm, ^ { - j + [—^— ) J , (2) 

where the minimization is over the locus of points j which defines the curve z(d, 0) toward 
each SN. A goodness of fit is computed by summing over all the SNe, 

X(P) 2 = EW)1 , (3) 

i 

with results shown at the top of Figure 2 for both the IRAS and ORS surveys. 

For the IRAS survey, we find (3iras — 0.40 ± 0.15 and for the ORS survey we find 
Pons — 0.30 ± 0.10. In both cases the value of x 2 at the minimum is within the expected 
tolerance, confirming gravitational instability as a valid model for the observed peculiar 
motions of SNe la. We find p°*^ s = 0.75 ± .38, in good agreement with the relative biasing 
of 0.7 derived from the correlation functions of ORS and IRAS galaxies (Fisher et al 1994). 

The predicted and observed peculiar velocity fields for the best values of j3 are shown 
in Figure 3 in the Local Group (LG) frame at the location of each of the 24 SNe la; 
these numbers are also listed in Table 1. Although the gravity predictions are derived 
independently from the observed peculiar velocities, the similarity of the two is remarkable. 
The leverage any single SN la measurement has in determining (5 depends on its location. 
Because the dominant feature of the flow is dipolar, SNe la along the axis of this dipole 
carry more weight than those whose radial motion is perpendicular to it. Although the 
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supernova data can be extended to greater depth, beyond 10,000 km s 1 , the sampling of 
the IRAS and optical surveys is inadequate to derive useful peculiar velocity predictions. 

3. Tests of Robustness 

To verify that our results are free from the vagaries of our supernova and galaxy 
samples, we performed bootstrap resampling tests on both data sets. This procedure tests 
the effects our choice of sample points and their uncertainties have on the estimates of f3 
by drawing new samples of data from our best estimate of the underlying population: our 
sample (Press et al. 1992). 

We first drew 200 sets of randomly chosen SNe la from the sample of 24 SNe la, drawing 
each object with a Poisson probability of expectation value of 1. We then subjected each 
set to our maximum likelihood estimator for (3, using the gravity field of the IRAS catalog. 
As seen on the bottom of Figure 2, the distribution for f3 preferred by the individual SNe 
la agrees well with the single likelihood distribution estimated from our sample. 

Similarly, we tested the robustness of the IRAS gravity field by generating 20 bootstrap 
resampled IRAS catalogs, as we did for the resampled SNe catalogs. For each resampled 
IRAS catalog, we generated full gravity fields for the range 0.1 < f3 < 1.2. The resulting (3 
values derived from the minimum \ 2 of the 24 SNe sample are also shown at the bottom of 
Figure 2. Again, this distribution is consistent with the constraints for f3 inferred from the 
X 2 of the original datasets. Thus we believe the estimates in §2 are reasonable. 



4. Discussion 
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Figure 3 demonstrates the remarkable consistency between the observed peculiar 
velocities of 24 SNe la and those predicted from the gravity fields of optical or IRAS 
galaxies using linear perturbation theory and the best value for f3. The excellent y 2 fits in 
Figure 2 confirm our simple model for the source of peculiar velocities while putting useful 
constraints on the mass density parameter, f3. 

The signals, seen in Figure 3, are largely dipole patterns revealing the motion of the 
LG relative to the SNe frame. They constrain the shear in the large scale velocity field 
induced by the gravity of galaxies within the sample's 10,000 km s _1 radius. The bulk of 
the CMB dipole signature does appear to have been generated within this radius. Velocity 
dipole patterns that match the gravitational dipole signature have been detected in several 
other surveys (e.g. Riess, Press, and Kirshner 1995b, da Costa et al 1996; Giovanelli et al 
1996), but flows inconsistent with the predicted gravity field for all (5 values have also been 
reported (Lauer and Postman, 1992, 1994; Davis, Nusser, and Willlick 1996). 

A more sophisticated, normal mode comparison of the IRAS gravity field to the 
velocity field derived from a sample of 2900 Mark-Ill Tully-Fisher galaxies within a limiting 
redshift of 6000 km s -1 shows inconsistencies for any value of (3 (Davis, Nusser, and Willick 
1996). This same procedure applied to the SFI catalog yields consistent gravity and velocity 
fields with Piras = 0.6 ±0.1. The POTENT reconstruction of the local density field from 
the Mark-Ill catalog recovers many observed features in the IRAS maps and estimates 
P = 0.89 ±0.12, but some inconsistencies persist (Sigad et al 1997, Dekel 1994). The failure 
to match the fields remains unexplained. 

By limiting the analysis to a redshift of 3000 km s -1 , Willick et al (1996) successfully 
applied the VELMOD algorithm to compare the IRAS gravity field to a sample of 838 
Tully-Fisher galaxies; a maximum likelihood analysis leads to (3 = 0.49 ± 0.07, is consistent 
with our results for the IRAS gravity field. Similar low values of the density parameter 
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emerge from the least action method applied to the flow field with 3000 km s _1 (Shaya, 
Peebles, and Tully 1995). These procedures are distinct and do not all suffer from the same 
biases. Thus it is encouraging that they are leading to consistent (and perhaps reliable) 
results. 

In Figure 3, the individual peculiar velocities in the SN map appear to be slightly 
larger than the velocities predicted at the same locations by the IRAS or ORS fields. We 
attribute this to measurement noise and velocity noise (small scale velocity flows) present 
in the SN data but not in the heavily smoothed IRAS or ORS velocity fields. Because of 
this noise, some SN velocity residuals do not match the sign (filling in Figure 3) of the 
IRAS predictions and/or the mode of the signs of the hemisphere in which they occur. 
But overall, the x 2 minimization for the supernovae takes these sources of error properly 
into account and produces the best matched f3 for the data sets. Our determinatiuons 
of (3 are not sensitive to our estimate of the velocity noise or SNe la errors, the latter 
limited by the small dispersion from smooth Hubble flow. As the distances for individual 
supernovae increase, the errors (in km s -1 ) increase, while for the IRAS and ORS data 
sets the smoothing becomes stronger and the amplitude of point to point variations gets 
smaller. Although more distant SNe carry less weight, we find the objects within as well as 
beyond 5000 km s -1 give consistently low values for (3. Nusser and Davis (1995) show how 
to avoid the effects of mismatched smoothing that appear here, but the present SN sample 
is too sparse to apply their methods. 

MLCS distances are precise enough to characterize the peculiar velocity field in the 
direction of each supernova. Yet, during this application we found instrinsic uncertainties 
still limit the precision of relative SNe la distances to no better than 5%. Future 
investigations into SNe observables may improve the precision of these distances or our 
understanding of their limitations. 



-10- 



As the supernova data set grows, the precision of this comparison between gravity and 
velocity will improve, and the methods of analysis that have been used on the galaxy data 
sets will become appropriate. No tool for mapping the peculiar velocity field has a brighter 
future. 
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Table 1: Peculiar Velocity Data 
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192.60 


51.40 
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-466. 


123. 


-407. 


-365. 


1996X 


310.20 


35.70 


1845. 


-96. 


149. 


-81. 


-204. 


1995D 


230.00 


39.67 


2000. 


-258. 
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-457. 


-191. 


1996Z 
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22.60 


2014. 


-427. 


289. 


-320. 


-71. 


1991M 
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2825. 


-381. 
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30.27 


3639. 
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187. 
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46. 
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226. 
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369. 
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327. 
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551. 
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-59. 
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564. 
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-29. 
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-49.00 
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252. 


221. 
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263. 


443. 


-156. 


62. 


1992ag 


312.50 


38.40 


7295. 


-489. 


984. 


-272. 


78. 


1992P 


295.62 


73.11 


7447. 


-519. 


539. 


-408. 


-108. 


1994Q 


99.60 


65.00 


8956. 


-790. 


648. 


-12. 


-632. 


1996C 


64.38 


39.68 


8872. 


-535. 


819. 


137. 


-477. 


1993ah 


25.90 


-76.80 
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-53. 
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183. 


277. 
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749. 


5. 
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180. 
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-260. 


23. 
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1000 2000 3000 

cz (km s _1 ) 

Fig. 1. — Velocity-distance relation in the vicinity of SNe la. This is an example of the 
velocity versus distance predictions of the IRAS gravity maps as a function of the mass 
density parameter, ft, in the direction of SNe la. The solid line shows the j3 — 1.0 prediction, 
with subsequent dotted lines showing the predictions for f3 in decreasing increments of 0.2. 
Overplotted is the measurement of the SNe la redshift, MLCS distance (in km s" 1 ) and 
distance error with a small scale velocity error of 200 km s _1 . Such plots contributes to the 
determination of (3 via equation (2). 
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Fig. 2. — Constraints on the mass density parameter, /3, and tests of robustness. The top 
panel shows the value of our \ 2 statistic as a function of j3, equation (3), reduced (divided) 
by the twenty-three degrees of freedom. Our best (3 is 0.4 from IRAS, and 0.3 from the 
ORS, with (3 > 0.7 and (3 < 0.15 ruled out at 95% confidence levels for the comparison to 
IRAS. Bootstrap resamplings (bottom panel) of the SNe la and IRAS galaxies, descibed in 
§3, validate our estimates of (5. 
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RAS (3 = 0.4 





Fig. 3. — Predicted and observed peculiar velocity fields in the Local Group rest frame. 
Filled/Open points represent SNe la with measured negative /positive peculiar velocities 
(bottom map). The top two maps show the peculiar velocities predicted by the gravity fields 
of the IRAS and ORS catalogs at the position of the SNe la for values of (3 which best fit 
the SNe la velocity field. Filled/Open crosses mark the direction towards which the Local 
Group is approaching/receding. The directions of increased/decreased cosmic microwave 
background temperature are indicated by +CMB/-CMB. 



